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Excitation of the odd parity quasi-normal modes of compact objects 
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The gravitational radiation generated by a particle in a close unbounded orbit around a neutron 
star is computed as a means to study the importance of the w modes of the neutron star. For 
simplicity, attention is restricted to odd parity ("axial") modes which do not couple to the neutron 
star's fluid modes. We find that for realistic neutron star models, particles in unbounded orbits only 
weakly excite the w modes; we conjecture that this is also the case for astrophysically interesting 
sources of neutron star perturbations. We also find that for cases in which there is significant 
excitation of quadrupole w modes, there is comparable excitation of higher multipole modes. 



0\ ■ I. INTRODUCTION AND OVERVIEW 

q\' 

Gravitational wave signals from perturbed black holes typically show waveforms with strong quasinormal (QN) 

ringing, damped oscillations at a period and damping rate characteristic of the mass and angular momentum of the 

hole. If the source of perturbations of the hole has no inherent time scale (as it would, for example, in the case 

fvj ■ of a particle scattered [|l|,|| or plunging into the hole 0) then QN ringing completely dominates the appearance 

of the waveform. Indeed, most perturbation computations show waveforms that contain only an initial transient 

C\| . and QN ringing. The meaning of QN ringing has been well studied in the mathematical context of perturbation 

^ ' calculations (Einstein's field equations linearized about a black hole background). QN ringing, however, has also 

CN found to dominate the waveforms of black hole processes computed with numerical solutions of the full nonlinear set 

of Einstein's equations ^^■ 

Unlike black holes, neutron stars have inherent time scales associated with the various possible modes of fluid 

^~s ■ motion. These fluid modes can be understood, and analyzed approximately with a semi-Newtonian approach: The 

^\ \ excitations and frequencies of the fluid modes are analyzed with Newtonian mechanics, and the fluid motions are 

On ■ treated as a source in the equations of weak field general relativity. The emission of gravitational waves, in these 

7j ' calculations, damps the periodic fluid oscillations of Newtonian theory, but due to the weak coupling to gravitational 

Q^i waves these fluid QN modes are weakly damped. By contrast, the QN modes of a Schwarzschild black hole are strongly 

' ■ damped, with damping times of the same order as their period. 

j^JT)! The sharp distinction between the nature of neutron star QN modes and black hole QN modes ended when studies 

J^ ' of neutron star models by Kojima |^ and Kokkotas and Schutz |Q revealed the existence of modes that had no 
Newtonian counterpart. These modes, called w modes (for wave modes) are general relativistic in their origin and 
are similar to black hole QN modes in that they are quickly damped. That w modes are distinct from fluid modes is 
^ ] particularly clear in the existence of odd-parity w modes. The quantities describing the perturbations of a spherical 
stellar model break up into two sets of quantities, called even and odd parity, which are not coupled by the perturbed 
field equations. The even parity set includes all the information about fluid motion S. The odd parity set includes 
no fluid perturbations. It describes purely relativistic modes that are perturbations only of the background geometry. 
Chandrasekhar and Ferrari JE ] showed that sufficiently compact homogeneous neutron star models had odd parity 
QN modes. Later studies fl^Jn] not only confirmed these results but also showed that these odd modes exist for all 
degrees of stellar compactness and that the less compact the star is, the more rapidly damped are its modes. Odd 
parity modes of neutron stars with realistic equations of state have also been computed recently [021 . 

In the next few years the vibrations of neutron stars may be of more than academic interest. Neutron star 
processes, in various forms, are among the most plausible sources of gravitational waves that may be within the reach 
of detectors like LIGO, VIRGO, GEO600, TAMA300 and others. Though the w modes are at frequencies above the 
range of high sensitivity for interferometric detectors ||l3| , they may be interesting sources for resonant bar detectors. 
In 1996, Andersson and Kokkotas [ [l4[ pointed out that much work on neutron star dynamics was still being done 
in Newtonian theory and these studies, in principle, excluded the fundamentally relativistic w modes. Andersson 
and Kokkotas particularly directed attention to the question that is most important about the w modes: to what 
extent are they excited in astrophysical processes. More specifically one should ask about the relative excitation of 
the semi-Newtonian fiuid modes and the fundamentally relativistic w modes in a real astrophysical process. 

The key difficulty here lies in the question of what is an acceptable model of a "real astrophysical process." Those 
processes which are of greatest interest as sources of detectable waves are the oscillations of a neutron star formed 
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in a supernova collapse or formed by the coalescence of two neutron stars as the last stage of binary inspiral. The 
numerical modelling of such processes with general relativistic codes is beyond near term computational capability, 
so "real" processes cannot be simulated. On the other hand it is not adequate to choose any convenient source of 
perturbations. 

Andersson and Kokkotas ( ||l4|), for instance, gave an example of the excitation of QN ringing of neutron stars 
by impinging gravitational waves. It is plausible that the relative excitation of w modes and fluid modes for such a 
source is not representative of how modes would be excited by astrophysical processes. This might be said also of 
any source of perturbations chosen arbitrarily. The model we investigate here avoids arbitrariness by using a specific 
physical process. We consider a nonrotating neutron star of mass M and radius R. An object of mass mo passes by 
the neutron star on a geodesic trajectory that is the general relativistic equivalent of a hyperbolic orbit. We make 
the "particle" approximation in which ttiq <C M and the field equations are expanded in the perturbation parameter 
mo/M. We study the dependence of excitation on the two parameters of the particle orbits, the particle energy 
and angular momentum. Neutron star models considered are limited to relativistic polytropes. In principle, our 
"incontrovertibly astrophysical" model would give a definitive answer to the question of the relative excitation of fiuid 
modes and w modes. But in practice the present paper gives only a partial answer and represents only a first step. 
Here, for simplicity, we consider only odd parity modes. The absence of coupling between odd parity perturbations 
and fiuid motions makes the odd parity case much simpler to analyze, but it also means that there is no excitation of 
fluid modes with which comparisons can be made. What can be achieved is a comparison of gravitational radiation in 
w modes, and gravitational radiation directly due to the moving particle. In this sense the present work also shows the 
parameter values for which w mode excitation is significant. These results are of interest by themselves and probably 
indicate the systematics of the even parity excitations. 

We present results for odd parity radiation from a range of particle trajectories and we use these results to speculate 
about the importance of w modes in astrophysical processes. We also demonstrate that the study of the excitation 
of strongly damped QN modes is better carried out with waveforms than with spectra. An additional point inherent 
in our results is that QN ringing of higher multipoles is excited as strongly as that of quadrupole modes, though this 
conclusion is probably a feature only of models involving perturbations of a neutron star by a relativistic particle. 

The remaining sections are organized as follows. Sec. O introduces the model problem. The wave equation descri bing 



the odd parity perturbations appear in Sec. [II, the analysis of the solution with Fourier transforms is given in Sec. IV , 
and the method of computing solutions i s pre sented in Sec.M Results are presented and discussed in Sec. VI, and 



conclusions are briefiy summarized in Sec. VII. Some relevant details of particle orbits in the Schwarzschild geometry 
are given in an appendix. Throughout the paper we use geometric units G ~ c— 1 and metric signature (+, — , — , — ). 
As this paper was in preparation we learned of two works dealing with excitation of QN modes of a neutron star: a 
paper by Tominaga et al. p5| on substantially the same problem investigated here, and a paper by Ferrari et al llfl] 
which considers the excitation of the even parity modes. Among the differences between these works and the present 
paper are an emphasis on waveforms rather than spectra as a mean to detecting excitation of QN modes, the study 
of high energy orbits and of odd modes of neutron stars of radius greater than 3M . 

II. THE MODEL: PERTURBATION OF A STATIC STAR BY A PARTICLE MOVING IN A 

SCATTERED GEODESIC 

In our model of gravitational interaction, we start with a static and spherically symmetric spacetime background 
metric 

ds^ = e''('')dt2 _ e^(r)dr'^ - r^[de^ + sin^ edip^] (1) 

describing both the interior and exterior of a star of ideal fluid of mass M and radius r ^ R. The stress energy of the 
fluid is 

Tap = (/9 + p)UaUi3 - pgaf3 (2) 

where p is the mass-energy density and p = p{p) is the pressure. The mass function m{r) is defined by e"^'-'"' = 
1 — 2m{r)lr, and the structure of the stellar interior is found by solving the hydrostatic equilibrium equations of 
general relativity 

dm 2 

— — — Airr p 
dr 

dv 2 dp 

dr p + p dr 



dp {p + p) ['m{r) + Airr^p] 

dr r[r — 2r7T,(r)] 



(3) 



We limit considerations here to stellar interiors that are relativistic polytropes p7[ . The polytropic equation of 
state has the simple algebraic form 



P/Po = (p/po)^+^/" , 



(4) 



A particular model is specified by fixing the polytropic index, n, (in general not an integer) and the central density 
and pressure po,po. With these fixed one can integrate to the surface (i.e, the radius at which p = 0) and find the 
value of the star's radius R, and mass M = m{R). The ratio M/R is an indication of "how relativistic" the stellar 
model is. A particularly simple class of polytropics is the n = constant density stars. For these models the metric 
coefficients are found to have the closed forms 



c'^ir) 



1- 



2Af 



1- 



2Afr2 



<R 



(5) 



-A(r) 



(l - 2MrVi?3 



<R 



(6) 



It is evident in (H) that 6"^^°) vanishes for M/R — 4/9. This well known limit to how relativistic a homogeneous model 
can be corresponds to infinite central pressure. Outside the star, the metric functions take their familiar Schwarzschild 
form. 



-A(r) ^ g^(r) ^ ^ _ (2Af/r) r>R, 



(7) 



A point particle of mass mo ^ M, starting from infinity, passes close to the star (but never enters it), is scattered 
and proceeds out to infinity. We treat the particle as a perturbation of the spacetime inside and outside the star, and 
we analyze the perturbations only to first order in mo/M. The radiation reaction is second order in this perturbation 
parameter, so we can ignore it and assume that the particle moves along a scattered geodesic of the unperturbed 
spacetime exterior to the star, i.e. of Schwarzschild spacetime. The metric of the perturbed spacetime can be written 
as 
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where the "(0)" index here and below denotes the background solution. The Einstein equations to first order in 
perturbations are 



SRap = 8Tr[STaf3 - -STg^ 



(0) 

Q/3 



ir(°)M . 



(9) 



where T^^'' = p — 3p. The perturbed stress energy has two contributions. One contribution is that of the fluid 
perturbations, a contribution we ignore for reasons given below. The other contribution is the stress energy of the 
scattered particle. We take the particle to be moving in the {9 — 7r/2) equatorial plane on a timelike geodesic in 
which its position as a function of coordinate time is given by z°'{t). With this notation the particle stress energy in 
the coordinates dl[) is given by 



^rj.c.p,parUcle ^ ,„^ ^ ^i^lj^ ^ [cos 0] ,5 [^7 - ^{t)]^'''^ '^^^ 
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Here R{t) and (f)(t) are respectively the Schwarzschild radial and azimuthal coordinates, as functions of Schwarzschild 
coordinate time, for the particle orbit. 



III. THE WAVE EQUATION GOVERNING ODD PARITY PERTURBATIONS 



A decomposition of perturbation tensors in tensor spherical harmonics (see Il8]) leads to a split of (ph into two 
sets of decoupled equations |]T9[ called even and odd parity. Fluid perturbations are purely even parity and couple 
only to even parity spacetime perturbations. Since we consider here only the odd parity (sometimes also called axial) 
perturbations, our analysis involves only the stress energy perturbations of the particle itself, as given in (10) 



We adopt the notation of Regge and Wheeler [gOJ and of Moncrief [gT| . In this notation the only nonvanishing odd 
parity metric perturbations, for a particular l,m multipole, are 

hte = -h'^''{r,t)id/sin0dip)Yi^ ht^ = h'o"'{r,t)ism0d/d9)Yi^ (11) 

hre^-h["\r,t){d/sm0dip)Yi^ K^ = h'{^{r,t)ism0d/de)Yi^ (12) 

hee = -h^^/sm'^0^ h'r{r,t){d'^ / sin 080 dip - cos 08/ sin^ 0d(p)Yim (13) 

hg^ = ^h^2"'{r,t){8^ / sm08(p8ip + COS08/80 - sm08'^ /8080)Yi^ . (14) 

For simplicity, for now we make the Regge- Wheeler [gO[ gauge choice, which for odd-parity means that we choose 
/12" = 0. By taking a combination of the odd parity perturbed vacuum Einstein equations, Regge and Wheeler |2C| ] 
showed that Qim (r, t) defined by 

Q,™(r,t)^e[''M-^Wl/2S!V^ (I5) 

r 

satisfies a sim ple homogeneous wave equation. The same combination of odd parity equations with stress energy 
perturbations |Lq| results in a wave equation with a source term, of the form 

%^-^ + {'a + l)-^+4..^[pM-PW]}^^'''^^-'5-(M). (16) 



Here r* is the usual tortoise coordinate 



dr 
which, in the vacuum outside the star, can be written as 



'^^ g[^(r)--A(r)]/2 ^ (--^7^ 



r + 2M log[r/2M - 1] + const . (18) 



The source term, Simir,t), is constructed from the stress energy ( [lOD of the particle. For a particle of mass mo, 
moving in the equatorial {0 = 7r/2) plane, with energy per mass E, and angular momentum per mass L, the l,m 
multipole of the source term has the form 

5..(M)^fl-^Vl-^)^-^^fl-^V + ^fl-^) (19) 

\ r J \ r J r ^ \ ^ J ^\ ^ J 

where 

Vim = DimS[r - Rit)] Cim = Cim5[r ~ Rit)] (20) 

with 

In the Schwarzschild spacetime outside the star dlQ) is the well known Regge- Wheeler equation |2fl] . 

IV. FOURIER TRANSFORM OF THE WAVE EQUATION 

The equation in ( |lq ) could be solved directly with a finite difference representation of the 1-1-1 partial differential 
equation. Our Cauchy data for this equation, however, correspond to the particle "starting" at infinity with no 
incoming radiation. These starting conditions are most simply handled with a Fourier representation of the problem, 
so we write Qim{r-,t) as a Fourier integral 

1 /"^ 
Qim{r,t) = — e-'^*Qim{r,uj)diu (23) 

27r ,_^ 



and transform ( Jig ) into a second order ODE for Qim{r,u)), 

^ '.''^'-HQi,n(r,uj) = -Sira{r,uj) (24) 
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where 



Sira{r,^)^ / Sira{r.ty'dt (25) 



is the Fourier transform of the source term. 

At infinity, the wave will be purely outgoing so that 

Qz„(r,c^) -> Ai™(w)e^-'-*, r* ^ oo . (26) 

It follows that at infinity, Qim is a function of retarded time u = t ~ r* given by 

Qirniu = t-r*) = ^ f Ai^{uj)e-'^^duj . (27) 

It is useful to note that the definition of Qim in terms of multipole components of the real quantities hij , and the 
property Yj*^ — (— 1)™Y;__,„ of the spherical harmonics, require 

AU^) = (-l)"M*_,„(-a.) . (28) 

The gravitational (odd parity part) power, radiated to infinity is given by p2| 

1=2 m=-l ^ ' 

Using ( |27|) and Parseval's theorem we obtain the odd parity part of the energy spectrum 

E §^l^.™(-)P . (30) 
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From the wavefunction Qimiu) we can construct the multipoles of the metric perturbations. Of particular interest are 
the transverse strains of gravitational waves. In an asymptotically flat gauge these are the perturbations hgg, h^^, hg^. 
To compute the relationship between Qim and perturbations in an asymptotically flat gauge it is useful to note that 
the quantity 
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has been shown by Moncricf pl| ] to be a gauge invariant combination of odd-parity metric perturbations in the 
Schwarzschild background. In the stellar exterior, in the Regge- Wheeler gauge, the perturbation Q^°'^^ agrees with 
our wave function Qim- It follows that Qim is identical to Qj^"'^ in the stellar exterior, and that in the exterior we 
can take ( |3l| ) to be the definition of ([1^) in an arbitrary gauge. Since the right hand side of (pi| ) is gauge invariant, 
we can evaluate the relationship between Qim and the metric perturbations in an asymptotically flat gauge, at large 
radius. In this case, the expression for Qhn in (pl|) reduces to 



1 r)/)'™ 1 . 

Qlnr = TT^ + 0{r-^) = --h2.lm + ©(r'^) . (32) 

2 or 2 

Here the dot over hij" indicates the derivative with respect to retarded time u. The last equality in (|32) follows from 
the fact that hif^ is a function of u aside from corrections higher order in 1/r. From ( p2| ) we have hi{^ and hence, 
with (y_l|) we have, for example, 

hee = h["'{r,t)(d'^ / sm9d9d(p - cos9d/ sin^ ed(p)Yim ■ (33) 



But this multipole component is complex. To get meaningful physical quantities we need to sum multipole components 
with ±771. In doing this it is useful to define: 
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For a given pair of multipole modes I, m and l^—m we can then construct real functions 
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= 4|m|r [Re (Q'") sin \m\ip + Im (Q'™) cos \m\ip] F[ 



^e" + ^e';™ = 4r- [-Re (Q'™) cos |m|</j + Im (Q'") sin |m|</7] F^ 



(34) 
(35) 
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V. COMPUTATIONAL IMPLEMENTATION 



A. Solution for Aim{i-o) 

If the background spacetime is due to a star, equation ( p4|) is to be solved for the boundary conditions of a 
regular solution at the center of the star and outgoing waves at spatial infinity: Qimif,^) -^ f''^^ for r — )■ and 
Qimi'f,'-^) — > e''^'' for r ~^ oo. The Green function solution is found in the usual way. We define 2/"^(r, w) and 
?/°"*(r, cj) as the homogeneous solutions of (p3) with asymptotic forms 



We then define the Wronskian of these two homogeneous solutions, an r independent constant, to be 



dr* 
With the above definitions, the Green function solution is 
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where the lower limit in the first integral, r^ = r*{rt), denotes the value of r* at the minimum radius (turning point) 
of the scattered orbit. In the limit of large r, a comparison of (p6|) with (140) gives us 
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With the notation for Cim,'Dim hi (II^) - ( pO[ ) the expression for Simir,uj) takes the form 
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Substituting this expression into (^11) , and using the properties of the S function, we get 



Almi^j) 



^{£.<e"'r'W.| 



2M \ Dim Clm 

^ R ) R^ R 



die 



.^tdf'^Dim 



dr R 



1- 



2M 



(43) 



For the geodesies eonsidered in this work (scattered geodesies, see appendix A), we can distinguish two phases of 
the motion: In the first, the particle comes from infinity and moves towards the star until reaching the turning point; 
it then returns to infinity in the second region. If we choose the origin of coordinate time and angle (p at the turning 
point, then for any time t — —t and angle ip — ~Cp in the first, inward, phase of the motion, there is a corresponding 
point t = t,ip = ip during the second, outward, phase. The symmetry of the inward and outward motions allows us 
to write the amplitude as an integral only over positive time. 
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where \Dim\ = -iDi,nexp{imip), \Cim\ = Cimexp{imtp) with Di^, Qm given b y (|2l| -P2D. Once the quantities Ai^i^) 



are computed, the waveform follows from (pq) and the energy spectrum from (pOf) 

If the particle moves outside a spherically symmetric black hole, not a star, then the amplitude can be obtained from 
(E3) with only two changes. The ingoing solution y"^^ (see p3|}, of the Regge- Wheeler equation must be substituted 
for y''°^ in the integral in ((Jif), and the Wronskian defined in (39h must be replaced by 
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B. Numerical method 



The first step in the determination of the amplitude and waveform is the computation of the regular homogeneous 
solution of the equation (|2j). The numerical problem can be divided in two parts: Integration of ( |24|) from the center 
of the star to its surface. And then from the surface to the point (outside the star) where we want to compute y'^'^^{r). 
Due to the presence of the centrifugal potential, which diverges at the center of the star, it is convenient to introduce, 
inside the star, the new function Zi, through 



yr^ = r'+'Zi{r) 



(46) 



and instead of 
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from a point tq, near the center of the star, using the starting conditions 
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until reaching the surface of the star, r — R where we switch back to y\'"^ and integrate 
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from r = R, with the starting conditions. 
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until the desired point r. 

Both (^) and (^) foUow from a rewriting of (|2J) in terms of the variable r. This avoids the numerically time 
consuming task of finding r{r*). 

If the unperturbed spacetime is a Schwarzschild black hole instead of a star, then instead of yj"^ we need to find 
yj"^, the purely ingoing (at the horizon) homogeneous solution of (24). This is done by substituting y[°® by yj"® in 



( pOf ) and integrating from a large negative value of r* with the "initial" condition j/J"^ ~ e *"'" . 

The computation of the Wronskian in (BS) or (Esl) is now straightforward. We simply compute y\°^ (or yj"^) at 
large values of r and use the asymptotic form of ?/; at large r 



yr(^-oo) 



1 



l{l + l) 



2iojr 



(53) 



to find the Wronskian by definition. 

One point worth of mention is that since 
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where F denotes the Gamma function, Aim 7^ 0, for a given Z, only when m = —I + 1, — ? + 3, ..., I — 1. Thus for I = 2 
only A21, A2-1 are non zero and for I = 3 only A30, A32, ^3-2 are non zero. 

A consequence is that if, as we will do in the next section, we restrict attention to the waveform seen by an observer 

in the plane where the particle moves, then hg~' = in (|3^) and the waveform will be polarized in the direction 

VI. NUMERICAL RESULTS AND DISCUSSION 

Our presentation of results starts with gravitational radiation in the quadrupole (I — 2) and octupole (/ = 3) modes, 
for a particle moving around a nonrotating black hole (see also |lj], ^). Though our main interest here is neutron 
star w modes, it will be useful to establish the black hole results for comparison. The case of a black hole will also 
be helpful in clarifying several issues relevant to the neutron star case. An additional motivation for the presentation 
of these results is to demonstrate that in cases in which I — 2 QN ringing is important, the I = 3 ringing tends to be 
comparable, both for black holes and neutron stars. 

It will be useful to understand the aspects of the emerging radiation that are associated with the motion of the 
particle, and those that may be ascribed to the influence of the background spacetime. (There is, of course, no 
fully satisfactory way of making such a distinction.) If the particle were in a periodic orbit with frequency Q, its 
gravitational perturbations would be functions oi ip — Qt, and /, m multipoles would have time dependence only at 
frequency mH.. For the scattering orbits one would expect that the characteristic peak of the spectrum would follow 
an approximate form of this rule with the particle's angular velocity dip/dt at the turning point rt taking the place of 
fl. In the discussion below we will use the expression "particle frequency" to denote this nonrelativistically expected 
peak frequency m dip/dt. Note that, as discussed above, an orbit in the equatorial plane is a source only for multipoles 
with m — I— 1,1 — 3, ■■■,—1 + 1. Thus, m will be 1 for quadrupole radiation and 2 for octupole. 

Figures |l| a and b show / = 2 and I — 3 spectra and waveforms for orbits with particle energy parameter E (energy 
per particle mass) of 2. For the particle motion in Fig. m the angular frequency dip/dt of the particle at the turning 
point is 0.18/Af so that the particle frequency is 0.18/Af for quadrupole radiation and 0.36/A/ for octupole. The 
waveform in Fig. nb shows more specifically that QN ringing is an important part of the waveform. 

The locations of the real parts of the least damped QN modes, uiM — .374 for I = 2 and ujM = .599 for I — 3, are 
indicated with arrows in Fig. n^a. The analogy with normal mode systems would suggest that the extent to which 
particle motion stimulates QN ringing increases the closer the characteristic particle frequency is to the QN frequency. 
In the case of Fig. Ha the (odd parity m = 2) particle frequency for I = 3 case is twice that for the (odd parity m = 1) 
I — 2 case, and the I — 3 QN frequency is around 60% higher than the / = 2 QN frequency. Both particle frequencies 
are roughly equal in the suitability for exciting the respective QN modes, and the plots show that QN ringing seems 
to be of roughly equal importance for I = 2 and / = 3. 



The presence of QN radiation is fairly clear in the waveform, shown in Fig. |l|b. This figure shows the gravitational 
radiation, a mixture of effects of particle radiation and QN ringing. It does not tell us what fraction of the radiated 
energy is in QN ringing, but that question probably cannot be given a clear formulation in any case |2J] . As an indicator 
of the importance of QN ringing, a plot of the waveform has a clear advantage, but it also has a disadvantage. To show 
the waveform Bq we must choose a direction in the 0, ^p 2-sphere. The curve in Fig. n^b corresponds to the equatorial 
plane {9 = vr/2) and to Lp = 7r/2. The appearance of the waveform in other directions will differ. As an indication of 
the complexity of the waveform dependence on direction, a quadrupole waveform in the equatorial plane is shown in 
Fig. ||. The radiation is that for a particle with E — 2 and L — lOM moving in the field of a constant density star 
with R = 3.125A/. The waveform amplitude, Re{Q2i) sin ip+ lin{Q2i) cos(p, is plotted above a cartesian plane with 
X = r* cos (/?, y = r* sin (p. The direction of the symmetry axis of the particle orbit is the straight line j/ = in the 
figure. Some symmetries of the radiation pattern can be seen in Fig. 0. Since the radiation is ttt, = 1 radiation, the 
amplitude at a given time (i.e. , at a given value of r*) must be antisymmetric under if —^ if + tt. At each radius in 
Fig. S, then, the variation in amplitude must have two nodes. It is interesting to note that the particle trajectory is a 
symmetric one (it has a line of bilateral symmetry) but the radiation pattern is not symmetric since there is a delay 
in time between the radiation given off by the initial and final portions of the particle trajectory. 

Despite this element of angular choice, the waveforms are the preferred tools for assessing the excitation of QN 
ringing in cases, like that of Fig. |l|, for which the QN mode "resonance" is broad. In the case of narrow resonances, such 
as the QN modes for ultracompact (and probably unrealistic) stellar models, the spectrum itself is the better indicator. 
This is clearly shown in Fig. ^ in which results are given for a constant density stellar model with R — 2.33M, very 
close to the lowest allowed limit R/M = 9/4. For this stellar model the four lowest damped / — 2 QN modes are 
u)iM = 0.281 + i0.00014, ivuM = 0.372 + i0.0048, uJmM = 0.458 + i0.025 and w/yM = 0.554 + i0.049. The results 
are given for / = 2 radiation from a particle with E ^ 2,L ^ lOM. Plots are presented for the waveform, in the 
direction 9 = 7r/2, if = 0, and for the energy spectrum. It is apparent in the latter that the first two QN modes of 
the star are excited and there are hints that all the first four might also be. The picture given by the waveform is 
rather more confusing due to the simultaneous presence of the two modes oj/ and loh with very small damping. The 
most obvious feature of the waveform, in fact is the beat frequency between these two almost periodic modes. For the 
case of this stellar model with several competing QN frequencies, and with sharp QN resonances, the nature of QN 
excitation is more clearly shown by the spectrum than the waveform. A second spectrum in Fig. ^ for a particle with 
a smaller energy of i? = 1.2 shows that the spectrum is equally useful for a particle orbit that is not highly relativistic. 
An interesting feature of the spectra is that the w modes appear as narrow spikes at the QN resonances. These spikes 
contain little energy, but they appear at the high frequency boundaries of broad spectral peaks very distinct from the 
particle radiation spectrum. 

The very sharp QN resonances appear to be unique to ultracompact stellar models. Since we will be concerned here, 
for the most part, with more-or-less realistic models, with broad QN resonances, we will rely on waveform pictures 
in our study of excitation of QN ringing. We will choose to clarify the presence of QN ringing in these waveforms 
by displaying them as linear-log plots of the magnitude of the waveform. In such plots QN ringing appears as a 
distinct sequence of bumps whose peaks lie on a straight line. The spacing of the bumps gives the real part of the 
QN frequency and the slope of the line gives the imaginary part. The distinctive nature of QN ringing in such a 
plot is shown in Fig. ]^ The waveform is shown in (a) and (b) respectively for the I = 2 and / — 3 multipoles of the 
orbit used for Fig. |l|. In (b) the waveform is also shown for an orbit with E = 3.7, L = 20M, and a turning point at 
rt = 3.767501Af . The pattern of bumps in part (b) has the same spacing and slope for both orbits, confirming the fact 
that the QN frequency is independent of the source. The higher oscillation rate and roughly equal damping for I = 3 
in comparison with / = 2 is evident in the patterns in (a) and (b). The QN frequencies found from these linear- log 
plots are ujM = 0.374 -|-i0. 090 for 1 = 2, and uiM = 0.600 -1-iO. 093 for I = 3. These are in excellent agreement with the 
known values cuM = 0.3736 -|- i0.089 for / = 2, and uM = 0.5994 + i0.0927 for I = 3. The combination of parameters 
{E = 3.7, L = 20M) for the comparison waveform in Fig. 0b was chosen because it gives excitation of QN ringing 
comparable to that ioi E = 2,L = lOM. It is worth noting that the particle frequency 2d(p/dt is almost identical for 
the two sets of parameters. 

In Fig. ra we present linear-log waveforms for a constant density stellar model with R/M = 2.5, an ultracompact 
stellar model that is astrophysically implausible, but that is useful as an illustration. The real part of the I = 2,m = 1 
waveforms are shown for particle orbits with the relatively low energy E = 1.2, and for a higher energy E = 2. For 
each energy one model is a low L orbit with a small turning point rt, and the other model has large L and large 
rt- In each case the large r* orbit gives far weaker QN ringing than does the small rt orbit. The long sequence of 
well defined bumps for the small rt models in Fig. g allow a graphical determination of the least damped I = 2 QN 
frequency tuM = 0.4123 -f- iO. 02199, to accuracy of better than 1%. 

Results for orbits around a constant density model with R/M = 3.125 are shown in Fig. |[ for two different 
values of L at each of two different values of E. For this model the least damped I = 2 QN mode has a frequency 
ujM = 0.455 + iO. 1386 with a real part slightly higher than for a black hole ( 0.39M^^) or an ultracompact R/M = 2.5 



model (0.41M~^). More important, the imaginary part of the QN frequency is large, around 50% larger than for a 
black hole and more than 6 times that for the ultracompact stellar model. The strong ringing part observed for the 
two smaller angular momenta can be approximated with precision better than 1% by the lowest damped QN mode. 

The patterns and trends in Figs. |g and |6| continue in Fig. [7|, in which waveforms are given for constant density 
models with R/M = 4, a value close to the typical radius of a neutron star. In this case the least damped I — 2 QN 
mode has frequency loM — 0.401 + i0.209, so the QN resonance is extremely broad. A consequence is that the pattern 
of QN bumps is less clear in the linear-log waveforms. For the low L, high particle frequency orbits, there is evidence 
of excitation of the lowest damped QN mode for both energies, but none for the high values of L. The energy spectra 
for these models are presented in Fig. H. For homogeneous stars with R/M > 4, the damping of QN modes will be 
even stronger than in Fig. 0, and QN ringing will be even less apparent in the waveform. 

The results presented in the previous figures are not sensitively dependent on the assumption of a homogeneous 
stellar model. Studies of polytropes with smaller stiffness show the same pattern of QN excitation as the n = Q 
constant density polytrope. A typical case is shown in Fig. Bl which displays the waveform when a particle with 
energies E — 2, A is scattered by a polytropic star of index n — 1 and with central density po = 5 x 10^^ g/cm'^. (In 
the equation of state p = KGp^, the parameter K is 10® m^). For this model the radius of the star is i? = 3.911Af, 
and the lowest damped I — 2 QN mode has the frequency uiM — 0.4603 + i0.1504. This mode is excited by both 
angular momenta considered in Fig. |9|b where the particle has a very relativistic energy, while it is only excited by the 
lowest angular momenta in Fig. 0a. From the graphs, the ringing frequency is computed to he ui = 0.458 + 10.1503, 
very close to the QN frequency. 

So far we discussed only excitation of I ^ 2 QN modes of stars, but as we have shown (see Figs. H and Ep), when 
the I = 2 QN mode of a black hole is excited by a scattered particle, the least damped I = 3 QN mode is excited 
with equal strength. To investigate whether this is also true for stellar models, we computed the octupole waveforms 
for a small set of homogeneous stars, with radii R — 2.5M, 3.125M and R = 4M. The results, displayed in Figs, [l^ 
and O, confirm that these compact stars have their I = 3 odd parity QN modes excited by scattered particles. The 
I — 3 results in these figures can be compared with corresponding I = 2 results. The orbits for Fig. |0p, are negligibly 
different from the orbits in Fig. |5|a. These figures describe radiation from orbits with the same energy and negligibly 
different angular momenta. The figures show that the I — 3 QN frequency, cuM — 0.598 + JO.0087, has a larger real 
part than that of the / = 2 ringing, but a slower damping time. The figures show that there is a similar range of 
angular momenta for which significant excitation of the quadrupole and octupole QN modes occurs. One can also 
compare Fig. ga with Fig. [lC|b. In both cases the particle moves outside an homogeneous star with radius R = 3.125M 
with the same energy, and again excitation is comparable. In this case, the I = 3 QN mode (wM = 0.691 + i0.143) 
has a damping time almost equal to the damping time of the I = 2 QN mode. A similar comparison can be made of 
Fig. [l| and Fig. ^ for stars with R = AM. 

VII. CONCLUSIONS 

QN ringing of w modes make a unique, high frequency, change in the appearance of waveforms. Even if QN ringing 
is only weakly excited it would in principle be easily recognizable and, as pointed out by Andersson and Kokkotas Q , 
would provide important information about the structure of the neutron star emitting the waves. Our results above 
show that / = 3, and higher multipoles, tend to be excited on the same order as / = 2 modes for scattering orbits. This 
is a somewhat unusual finding since gravitational processes tend to be dominated by quadrupole radiation. Significant 
higher multipole radiation may be specific to scattering by relativistic orbits of a "particle" (with negligible angular 
extent), and may not be generalizable to more astrophysically plausible sources. If these higher I modes do turn out to 
be as highly excited as the I — 2 mode, the modes will all be mixed together in the waveform recorded by a detector; 
a detector at a single location cannot "project out" individual multipole moments as we have done in the figures 
presented above. The signal from a neutron star could then be a highly complex pattern that in principle contains 
a great deal of information about the neutron star. In practice, the signal to noise ratio, and the breadth of the QN 
"lines" would make it very unlikely that QN mode spectroscopy could reveal the details of neutron star structure. 
Another practical consideration is that the w mode frequencies (typically 10^ Hz) are well above the upper limit 
(several hundred Hz) of laser interferometric detectors, the detectors that are most capable of detecting waveforms. 
The w mode frequencies, of course, scale inversely as the mass of the neutron star. But this is of little consequence, 
since neutron star masses are constrained to a very small range. By contrast there is at least the possibility for 
black holes of large enough mass (hundreds of solar masses) that QN ringing will be within the optimal bandwidth of 
interferometric detectors. 

From the representative waveforms and spectra presented we can infer that the lowest frequencies of gravitational 
radiation are negligibly affected unless QN ringing is very strong, and hence for a laser interferometric detector, an 
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analysis omitting general relativistic effects would be satisfactory. On the other hand, for a detector with sensitivity 
extending to hundreds of kilohertz, the w mode contribution could be very important. The ringing of w modes 
provides a high frequency component of the spectrum that could contain more energy than that due "directly" to the 
particle. For that reason it is of some interest to look through the neutron star results to make inferences about the 
conditions necessary for w modes to be strongly excited. It is tempting to conclude that the excitation of w modes 
is largest when the particle frequency is closest to the real part of the QN frequency. This conclusion, however, is 
somewhat uncertain. The particle frequencies are always below the QN frequency, so a higher particle frequency is 
the same as a particle frequency that is closer to the QN resonance. The most important determinant of the particle 
frequency is the turning point radius r^, and smaller rt will produce higher particle frequency, and hence a particle 
frequency closer to the QN resonance. But a smaller rj also means that the source of the disturbance to the neutron 
star is closer and hence is better able to excite QN ringing. 

The difference between these two interpretations is important if we try to generalize from the above results to 
neutron star processes that arise more naturally astrophysically. In this connection it is helpful to turn attention 
to the (astrophysically implausible) ultracompact stellar model of Fig. ^ . In the R/M — 2.33 model the two least 
damped QN modes have almost negligible imaginary parts, so the resonance picture would imply that excitation 
could be stimulated only by the high frequency tail of the particle radiation, which is negligible at the QN "resonant" 
frequency (the real part of the QN frequency). In spite of this the spectrum in Fig. ||b is far from a particle spectrum. 
There is in fact more radiation in the high frequency part of the spectrum than in the low. This implies that care 
must be used in applying the ideas of normal mode systems to QN modes, and that excitation may be much larger 
than a simple resonance model would suggest. 

For the models that are astrophysically most relevant, those with R/M > 4, the model results suggest that even 
the source most tightly coupled to the QN modes (i.e., even the smallest turning point) can stimulate only moderate 
ringing. This may be associated with the very short damping times for these stellar models and may not apply 
to stellar models with longer damping times |^,^. If we disregard such possibilities, and assume that the QN 
frequencies of stellar models are similar to those of polytropic equations of state, we infer that w mode effects on 
the radiation from scattering orbits are not of crucial significance. We conjecture that supernova core collapse or 
the coalescence at the end point of binary inspiral will couple no better to w modes than do scattering orbits. For 
R/M > 4 models, then, w modes will be of crucial importance only if the source of spacetime perturbations does not 
stimulate any radiation due to fluid motion, as is the case for excitation of ringing by an odd parity gravitational 
wave. 
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APPENDIX A: TIMELIKE SCATTERED GEODESICS OF SCHWARZSCHILD SPACETIME 

In this appendix we present some of the properties of the family of geodesies of Schwarzschild spacetime considered 
in this work. Additional details can be found in [ p7[ . 

A particle with a given energy E and angular momentum L (both per unit rest mass), moving in one of these 
geodesies starts its motion at radial infinity with initial velocity Wqo = y/E'^ — 1/E, reaches the turning point rt and 
returns to infinity. 

Scattered geodesies require E >= 1. The higher the energy of the particle the more relativistic it is (when E = 1 
the particle falls with zero initial velocity) . For a given energy there is a minimum angular momentum if the particle 
moves in one of these geodesies (see ex. M). It is 



ii E > 1 and 



L^ 27E^ - 36E^ + 8 + .y{27E^ - 36E^ + Sf + 64:{E^ - 1) 

M2 > 2(^2 _ 1) ( ) 



Jf>' (^2) 
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if £■ = 1. If the particle has smaller angular momentum and still moves along a geodesic than it will plunge into the 
star (or black hole). The equality happens for a particle moving in a circular geodesic. 

The turning point is the largest of the two positive roots of the cubic equation (the third root is always negative) 

r^{l - E^) - 2Mr^ + rL^ - 2ML^ = (A3) 

We can find it numerically by noting that it is always located in the region 



L2 / / 12M2 \ L^ + y/L^ - IMPL^ 



1 - \/l - ^^1 <rt< J77 (A4) 



and that the other positive root is always located outside this region. For a given energy, the smaller the angular 
momentum of the particle is, the smaller r^ is. However the turning point is always 

n > 3M (A5) 



In the special case E — 1, (A_3) becomes a quadratic equation and the turning point can be written in analytical form 



L^ + VL^ - 16M2L2 
'^ ^ 4M (^') 

Clearly rt > 4:M in this case. 

Once we know the turning point for a given energy and angular momentum of the particle, we can determine the 



other two roots of (A3) easily. We need first to compute the parameters e,p (see for ex. p8[). 



_ L^rt ~ 2Mr1 + Ly/^UL'^M'^ + AL'^Mn + LV? - 16M2r2 
"" ' 2Af(L2 + ^2) (A7) 

P^^ie + 1) (A8) 



(note that e = 1 if £' = 1) after what we can rewrite (A3) like 

{r ~ ra)ir - rsec){r - rt) = (A9) 

with the three roots given by 

^0 — e-l' ^sec — p_4 , rt — ^^^ l^iUj 

If we choose the origin of coordinate time T, angle ip and proper time r at the turning point then the value of these 
quantities at a given value of the radial coordinate r will be determined by the geodesic equations, 

T{r) = ± ^ ^ / , (All) 

VE^^TJr, (1 - 2Af/r)v/(r - ro)(r - r,ec)(r - n) 

T C°° rlr 

^(r) = + / (A12) 

VE^ - 1 Jrt y^r{r - ro){r - rsec)ir - n) 

T(r)=±^^=/ , (A13) 

with the — sign holding for the part of the motion from infinity to rt and the + sign for the part of the motion from 
rt to infinity. 

Although the above integrals are finite, their integrands diverge at the turning point, making difficult their numerical 
computation. A possible remedy is to make a change of variable that makes the integrands finite over all the range 
of integration. For example 

z = ^/r-rt (A14) 

leads to 
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for ii^ > 1 and to 



Tir) = ± , / , (A15) 

2 /•2('-) r^/'^dz 

V-E-^ - 1 Jo \/[r -rQ){r-rs^c) 



^/'^dz 

(A18) 




M 7o (1 - 2M/r)Vr-rsec 

V A7 Jo Vr(r - Tsec) 

nr l-z(r) 3/2J 



if S= 1. 
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FIG. 1. Gravitational radiation from a particle scattered by a nonrotating black hole. The particle's energy per mass is 
E = 2 and its angular momentum per mass is L = 10, so that the orbital turning point (minimum radius) is rt = lOM/3. 
(a) The energy spectrum for I = 2 and I = 3 radiation. Arrows indicate the oscillation frequency (i.e. , real part of the QN 
frequency) for the least damped 1 = 2 mode {uAI = .374) and the least damped 1 = 3 {ujM = .599) modes, (b) The 1 = 2 
waveform in the equatorial {6 = n/2) plane. The real part of the wavefunction Q21 is shown; this corresponds to the waveform 
seen at azimuthal angle tp = n/2. 
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FIG. 2. The full waveform as seen by observers in the plane were the particle moves {6 = 7r/2). The background spacetime 
is a homogeneous star of radius R = 3.125M and the energy of the particle is E = 2 with an angular momentum L = lOM. 
The waveform as seen in the same plane in the direction ip = 7r/2 is shown below in Fig. 5(b). The middle part (spiral part) 
represents the ringing part of the wave. The inner (center) the tail. The outer part is the initial burst. In the figure, the (p = 
symmetry axis of the particle orbit is the line y = r* sin ip = 0. 
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FIG. 3. Gravitational radiation from a particle scattered by a homogeneous star of radius R = 2.33M. (a) Absolute value 
of the imaginary part of Q21 in a linear-log scale for a particle with parameters E = 2,L = lOAI and turning point lOM/3. 
(b) The energy spectra for the same particle and for a particle with a lower energy {E = 1.2, L = 6A1 and turning point 
n = 5.163M). The four lowest damped / = 2 QN modes, with frequencies ujiM = 0.281 + iO. 00014, cunAd = 0.372 + i0.0048, 
iUiiiM = 0.458 + i0.025 and ujivM = 0.554 + i0.049 are indicated by arrows. 
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FIG. 4. In (a) the absolute value of the I — 2 waveform of figure l.b is plotted in a linear-log display. The / = 3 waveform 
in the plane 9 = 7r/2 and direction i^j = 0, is shown in (b) along with the / = 3 waveform for an orbit with E — 3.7, L — 20M, 
and rt = 3.768M. 
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FIG. 5. The real part of the I = 2,171 — 1 odd parity waveform Q21, for a constant density star with R = 2.5M. The 
waveforms in (a) are for two orbits with E — 1.2. The low L orbit has L = 5.5M, rt = 4.047M, and a particle frequency 
dif/dt = 0.142M"i; the high L orbit has L = lOM, n = 11.633M, and a particle frequency dip/dt = 0.051M"^ The waveforms 
in (b) are for two orbits with E — 2. The low L orbit has L = lOM, rt = lOM/3, and a particle frequency dip/dt = 0.180M~ ; 
the high L orbit has L — 30Af, rt — 15.862M, and a particle frequency difi/dt = 0.052Af~^. The least damped {I — 2) odd 
parity QN mode of this star has frequency ujM = 0.412 -I- i0.0219. 
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FIG. 6. The real part of the I = 2,m = 1 odd parity waveform Q21, for a constant density star with R — 3.125M. The 
waveforms in (a) are for two orbits with E = 2. The waveform with pronounced QN ringing has L = lOAI, rt — lOM/3 and a 
particle frequency difi/dt = O.ISOM"^ ; the waveform with weak QN ringing has L = 25M, rt = 12.94M and d(p/dt = 0.063M"\ 
The waveforms in (b) are for particle orbits with E = 3. The waveform with strong QN ringing corresponds to a particle orbit 
with L = 16A/, rt = 3.729M and d(p/dt = 0.178M~'^; the waveform with weak ringing corresponds to a particle orbit with 
L = 45M, n — 14.66A4" and dip/dt = 0.0603M~^. The strong ringing part can be fitted with precision better than 1% by an 
oscillating damped function with the frequency of the lowest 1 = 2 QN mode which is ujM = 0.455 + i0.1386. 
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FIG. 7. The real part of the / = 2, 7n = 1 odd parity waveform Q21, for a constant density star with R = 4Af. The 
waveforms in (a) are for two orbits with E = 2. The waveform with pronounced QN ringing has L = IIM, rt = 4.34M 
and a particle frequency dip/dt = 0.157M~^; the waveform with no evidence of QN ringing has L = 25A/, rt = 12.94M 
and dip/dt = 0.0631i\f ~ . The waveforms in (b) are for particle orbits with _B = 4. The waveform with strong QN ringing 
corresponds to a particle orbit with L = 23M, rt = 4.259M and dip/dt = 0.1682A/~^; the waveform with weak ringing 
corresponds to a particle orbit with L = 35M, rt = 7.71M and dip/dt = 0.109Af~^. For a constant density star with R = A 
the least damped odd parity [I = 2) QN mode has frequency u)M = 0.401 + JO.209. 



10-' 


r 


\ 


\ 


\ 


•J 


r 


- 


■... 




E=2,L/M=11 

E=2,L/M=25 


10"' 










\ 


\ 


\ 


\ 




V^^^^ " 


10-' 
















\ 
\ 




^^^ 


10-' 
















\ 


\ 
\ 


"n 




FIG. 8. The energy spectra for the case depicted in the previous figure. The background spacetime is a homogeneous star 
of radius R = AM. 
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FIG. 9. The absolute value of the imaginary part of the odd parity quadrupole waveform Q2\ is shown for a polytropic star 
of index n = 1 and K = IQQKm? . The central density of this star, po = 5 x IQ^^g/cm? and the radius is _R = 3.91Af. In 
(a) the scattered particle has energy E = 2 and two angular momenta, L — IIM (turning point 4.344M,d<^/di — 0.157M~'^) 
and L = 25M (turning point 12. 94SM, dip /dt = 0.063M~^) are considered. In (b), a higher energy {E = 4) is considered and 
also the waveform for two angular momenta are shown: L = 23Af (turning point 4:.259M,d(f/dt = 0.168M~^) and L = 35M 
(turning point is 7.71M,dip/dt = 0.109M"^). 
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FIG. 10. The absolute value of the real part of the odd parity octupole waveform Q32 is shown for a homogeneous star 
of radius R = 2.5M in (a) and of radius R = 3.125M in (b). In (a) we consider a particle with small energy {E = 1.2) 
and two angular momenta: L = 6M (turning point is 5.163M,2d(p/dt = 0.23Af ~^) and L = lOAf (turning point is 11.633M, 
2d(p/dt = 0.102Af~^). The amplitude of the wave for the far orbit is so small that the ringing observed (at about twice the 
frequency of the QN mode) is probably numerical error. In (b) we consider a particle with energy {E — 2) and two angular 
momenta: L = lOM (turning point is WM/3M,2dip/dt = 0.36M"^) and L = 20M (turning point is WM,2d(p/dt = 0.16M~^) 
moving outside a star with radius 7? = 3.125M. There is, for both angular momenta, excitation of a QN mode of the star, 
that most likely is the lowest damped I = 3 odd mode (from the graphic for the smaller angular momentum orbit we get 
L^A// = 0.691 +iO. 143). 
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FIG. 11. The absolute value of the real part of the odd parity octupole waveform Qzi. is shown for an homogeneous star 
of radius K = 4M. In (a) we consider two sets of parameters for the particle: E = 2,L — IIM (turning point is 4.343M, 
2d(p/dt = 0.315M~^) and E ^ 3,L = 17M (turning point at 4.25M, 2d(p/dt = O.SSSM"^). In both cases ringing at the 
same frequency is visible (that frequency is approximately cu — 0.64 + iO.23. In (b) the particle has energy E = 4 and we 
consider two angular momenta: L — 23 (turning point rt — 4.259M,2dip/dt = 0.336M^^) and L = 35M (turning point 7.71M, 
2dip/dt — 0.22M~^). In both cases there is ringing. The oscillation is clear enough for the smaller turning point to allow us to 
obtain its complex frequency: uiM = 0.64 + iO.23. 



18 



